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Abstract. We explore the super-horizon evolution of the two-point and three-point cor- 
relation functions of the primordial density perturbation in randomly-generated multi-field 
potentials. We use the Transport method to evolve perturbations and give full evolutionary 

O histories for observables. Identifying the separate universe assumption as being analogous 

to a geometrical description of light rays, we give an expression for the width of the bundle, 

Ci thereby allowing us to monitor evolution towards the adiabatic limit, as well as providing a 

useful means of understanding the behaviour in /nl- Finally, viewing our random potential 

I as a toy model of inflation in the string landscape, we build distributions for observables by 

evolving trajectories for a large number of realisations of the potential and comment on the 

prospects for testing such models. We find the distributions for observables to be insensitive 

to the number of fields over the range 2 to 6, but that these distributions are highly sensitive 

to the scale of features in the potential. Most sensitive to the scale of features is the spectral 

index, with more than an order of magnitude increase in the dispersion of predictions over the 

range of feature scales investigated. Least sensitive was the non-Gaussianity parameter /nLj 

7—i which was consistently small; we found no examples of realisations whose non-Gaussianity is 

^ capable of being observed by any planned experiment. 

• ^ 
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1 Why random potentials? 
1.1 First reason 

With the prospect of improved data from the Planck mission fast approaching, there has been 
a lot of interest in finding infiationary models exhibiting specific observable footprints. Large 
non-Gaussianity signals peaking at various shapes is one example, another being features in 
the primordial power spectrum. While this is a crucial step towards understanding what 
observables are specific to a particular model, often the set-up can be somewhat contrived 
and to gain understanding as to whether such behaviour is a general feature of the model, 
one may need to invoke Monte Carlo techniques. In such a situation it may be helpful to 
employ some element of 'randomness' at the level of the construction of the model. 

For example, a popular model of inflation coming from string theory is Dirac-Born— 
Infeld inflation. The DBI Lagrangian is 



c = -r(0) Ji + r(0) - v{^), (1.1) 
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where X = —^g^yd^^(j)d^(j) and T{(j)) is the brane tension. If the D3-brane velocity approaches 
its hmiting speed 

IX 



then a period of inflation can occur. This model has interesting observational consequences 
as the sound speed can become small and hence, since ~ l/Cg, lead to large equilateral 
. However the above Lagrangian also admits inflation by other means. In Ref. [1] a rather 
sophisticated model of brane inflation was investigated, where to simulate the effect of the 
bulk in different compactifications, random coefficients were used. In this set-up, conditions 
for DBI inflation were never encountered; instead inflection-point slow-roll inflation was vastly 
more common. We therefore see that while the DBI effect certainly gives an interesting 
observational footprint, there is no reason to believe this is a generic feature of brane models 
of inflation.^ 



1.2 Second reason 

On a more ambitious note, string theory seems to predict the existence of a landscape [2, 3], 
where, in the low-energy approximation, different regions may be characterised by the val- 
ues of a large number of scalar fields. The consequence of this is that we have some very 
complicated potential V{(j)i, ■■■,4>d) with a large number of minima each corresponding to a 
different metastable vacuum energy. This implies that instead of trying to predict the values 
of observables, we should be trying to predict probability distributions for them. Indeed, as 
we will now discuss, this is the case not just for string theory but for any model with multiple 
light fields. 

Most work on the consequences of a landscape has focussed on the measure problem 
(see Ref. [4] for a recent overview) but if the observational consequences of such a model 
are ever to be understood, then there are other challenges to contend with. In order for a 
landscape model (any model where the scalar potential has more than one minimum, or for 
the purposes of this discussion, even just one minimum but multiple fields) to be predictive, 
three questions need to be addressed: 

1. What are the statistical properties of the landscape 

2. What are the selection effects from cosmological dynamics 

3. What are the anthropic selection effects 

The measure problem relates to the question of handling the numerous infinities which 
turn up. Taking the example of slow-roll inflation,^ any model of multi-field inflation suffers 
from an uncountably infinite set of choices for initial conditions. In general one needs to 
assume that, one way or another, at some point our region of spacetime experienced field 
values displaced from our local minimum. This corresponds to a single realisation of initial 
conditions (plus quantum scatter) that gave an anthropically suitable inflationary trajectory, 
which subsequently found its way to our local minimum. In order to make predictions we 

^As pointed out in Ref. [1], this result is not conclusive since, rather importantly, their investigation did 
not go all the way to the tip of the throat. Nevertheless we feel this example illustrates the point in hand. 

^Most discussion in this area focuses on the scenario of inflation coming from tunnelling between metastable 
vacua but as we discuss here, the problem is much more general than that, affecting even the most pedestrian 
of inflationary set-ups. 
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need to ask what proportion of the whole universe finds itself in this situation, i.e. what 
proportion of an infinite space finds itself in one of an infinite set of initial conditions. The 
ratio is ill-defined without a measure. 

However overcoming this formidable task is not the end of it. Even with a solution to 
this measure problem we are still left with a considerable challenge. A solution to this issue 
is likely not to give us a specific set of initial conditions for a given model but a probability 
distribution for them. If all we can hope for is a statistical description of initial conditions, 
then in turn we only have a statistical description of infiationary trajectories and so, rather 
than calculating single values for observables, we should be calculating their distributions! 
The shape of these distributions will in part be determined by the model. This last point 
can, at least in some respects, be studied in its own right without a detailed knowledge of the 
string landscape or the measure problem. In this paper we take inspiration from the string 
landscape and study characteristics of these distributions in the context of a potential with 
multiple fields, containing a large number of vacua. 

An early study of the possible consequences of this landscape picture for slow-roll 
inflation was carried out by Tegmark [5], who generated a large number of random one- 
dimensional potentials and explored the inflationary outcomes. In Ref. [6] we extended this 
to two fields to investigate the effect of entropy modes on super-horizon evolution. As al- 
ready mentioned, in Ref. [1] a similar analysis was done for a six-field brane inflation model 
with random terms arising in the contribution coming from the bulk, where although the 
super-horizon effects were not analysed, both reassuringly and rather excitingly, qualita- 
tively similar emergent behaviour was identified to that found in Ref. [6]. In this paper we 
further extend our work in Ref. [6] to a larger number of fields and a broader range of po- 
tentials, as well as obtaining results for the non-Gaussianity /nl- The aim is to gain insight 
into the origin and limits of emergent behaviour. 

2 Architecture 

2.1 Not the most general Fourier series 

We construct our potential following an approach similar to Refs. [5, 6]. We use a random 
function of the form 



where (j) is the vector (pi with i running from 1 to d, as is the vector k and and mh are 
the vertical and horizontal masses respectively. The summation in ki means d summations 
take place where in each case the summation runs from ki — 1 to ki — ^max* 

The amplitudes 

= aki^...,kd independent Gaussian random variables with zero mean and standard 
deviation 



Due to computational limitations we cannot make /cmax sufficiently large to see the effect of 
the central limit theorem come into play. So rather than having a variance of order unity, we 
control the variance in the above manner. While this is a helpful thing to do computationally, 
one needs to bear this in mind when considering principles of effective field theory. We will 
return to this discussion shortly. 




(2.1) 



a = e 



(2.2) 
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Figure 1. Example trajectories for our potential (blue) and Fourier series (red) for two-field model 

with femax = 3. 



The potentials we simulate are periodic with periodicity scale 27rmh, and we can only 
expect reasonable results if the field trajectory spans a distance in field space less than the 
periodicity of the random function. This turns out to always be the case. 

Note that summing the potential this way means we are not using the most general 
Fourier series. As shown in Fig. 1, by restricting the summations over each fcj to non- 
negative values we sacrifice statistical d-spherical symmetry but in doing so we are able to 
build an observationally indistinguishable potential out of a fraction of the number of terms 
(see Table 1 for the number of terms in the series for various d and A^max values). 

2.2 Experiment-specific considerations 

As discussed, any model of inflation where the potential has multiple minima predicts a prob- 
ability distribution for the cosmological parameters. We wish to compute this distribution 
for various potentials of the above form. To do this we perform the same experiment as that 
performed in Refs. [5, 6]: 

1. Generate a random potential T^(</>) and start at (/> = (0,0). 

2. If V{0, 0) < then reject model, otherwise evolve to find the field trajectory. 

3. If model gets stuck in eternal inflation, reject. 

4. Once the model stops inflating, if the number of e-folds of inflation < 60 we reject 
as insufficient inflation occurred, otherwise calculate observables. 

5. Repeat steps 1-4 many times to obtain a statistical sample. 

6. (Change some assumptions and do it all again.) 
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Our Potential 


Fourier Series 


Fields d 


2 


2 


4 


4 


6 


6 


2 


2 


4 


4 


6 


6 


kmsux 


3 


5 


3 


5 


3 


5 


3 


5 


3 


5 


3 


5 


# of terms 


18 


50 


162 


1250 


1458 


31250 


98 


242 


4802 


29282 


235298 


many 



Table 1. Summary of how the number of terms in the potential changes with the truncation k, 
and number of fields d for Fourier series potential and our reduced version. 



Note that due to our potential being statistically invariant under translation, generating 
multiple realisations of the potential and starting at the origin is equivalent to taking a 
single realisation and scanning over initial conditions. 

In Ref. [6], taking the final minimum as the ultimate vacuum energy, to give an ap- 
proximately anthropically suitable solution [7] we had an additional cut stipulating the final 
vacuum energy must be positive to avoid subsequent collapse. This, in conjunction with the 
rejection of eternally inflating vacua, was found to be an extremely severe cut, in some cases 
reducing the proportion of otherwise viable solutions from 0.06 to more like 2 x 10~^. In this 
paper we abandon this cut, to enable us to explore more featured potentials which would 
not otherwise be computationally accessible. We found this to be of little consequence for 
observables. An explanation for this is that the two models may differ only in the nature 
of the post-inflationary evolution of the trajectory, which has no effect on the evolution of 
observable quantities. 

The other consideration regarding the experimental set-up is at what value to set the 
vertical and horizontal mass scales and m^. The vertical mass has little dynamical impact 
and only affects the amplitude of the observed power spectrum by a factor and not other 
observables. For this reason, rather than fixing we adjust it on a case-by-case basis such 
that the amplitude at horizon exit is = 2 x 10~^. 

The horizontal mass is more interesting. As previously discussed, motivated by the 
aim of minimising the number of terms in the potential for a given dynamical behaviour, the 
random coefficients are chosen in such a way as to make the potential essentially insensitive 
to truncation. This set-up means that m-h is our key parameter in adjusting how featured the 
potential is. As the examples in Fig. 2 show, adjusting the scale of features affects the length 
scale A(p of the inflationary distance in field space. We will discuss motivation from theory 
for this length scale next, but it is important we understand its implications for predictability 
and thus we shall be showing results for a range of m-h values. 

2.3 Basic considerations from effective field theory 

Generally one expects that the inflationary potential can be well described by an effective the- 
ory containing non-renormalizable contributions coming from integrating out massive fields.'^ 
For instance, for a single-field model one can write 

V{ct>) = Vo + -m'ct? + M,/.3 + -Ac/.^ + A^M^i j , (2.3) 

where the terms in the summation are non-renormalizable. One expects the masses in the 
summation to be at or even well below the Planck mass as, in analogy to the argument 

^Discussion along these lines can be found in Ref. [8] and of course Ref. [9] but we would particularly like 
to thank Liam McAllister and Sam Rogerson for very helpful clarifications and additional comments on this 
matter. 
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Figure 2. Example of two-field potentials with nih = i5.8Mpi and rrih = 2.0Afpi respectively. 

from W~W scattering for a Higgs around 1 TeV, tliere needs to be sometfiing to unitarize 
graviton-graviton scattering. Tliere is no good reason to assume the inflaton does not couple 
to these extra degrees of freedom. To do so is to make a strong assumption about quantum 
gravity which is hard to justify, and thus we expect ~ 1. If we categorise inflation models 
as large field, \A(f)\ » Mp\, medium field |A(/)| ~ Mpi and small field |A0| <C Mp\, then 
this sort of reasoning indicates small and perhaps medium field models should be considered 
more realistic as terms in the summation are suppressed, while to have a large-field model, 
one needs to justify additional symmetries to protect the flatness of the potential against 
the otherwise increasingly large series contributions. Crudely speaking we can think of our 
choice of rrih as corresponding to a decision on what energy scales we are integrating out. 

Finally, we would like to consider the number of fields to be included in the model. 
Historically a lot of focus has been given to single-field models simply because they are the 
most basic inflationary set-up, but this is not what is best motivated from the field theory 
perspective. As already mentioned, a single-field model occurs when one degree of freedom 
is much lighter than all the others. This means one can integrate out the other degrees of 
freedom provided they are sufficiently massive, but there is no good reason to believe this 
is necessarily the case. For example, in string theory the contributing massive fields include 
stabilised moduli. Work on flux compactifications is still very much in development but 
typically masses correspond to around the Hubble scale. This strongly motivates models 
with tens if not hundreds of active fields [10]. 

It therefore seems quite reasonable to model the final inflationary phase in a landscape 
as a truncated d- field Fourier series with random coefficients, provided we are dealing with 
small- to medium-field models. However, for computational reasons we are forced to work 
with something less realistic. Ideally we would work with more fields, and push to smaller 
field excursions than we will be working with. For the purposes of our investigation we will 
at times be working with inflationary trajectories that not particularly well motivated as 
genuine models of inflation, yet we still find them to be quite informative when it comes to 
understanding inflationary dynamics. 
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3 The Transport equations 



In this paper, we improve on our previous work [6] by calculating the perturbations using the 
Transport method of Mulryne et al. [11, 12]. This gives improved computational efficiency 
for the power spectrum, while still including all isocurvature effects, and additionally allows 
us to compute the non-Gaussianity parameter /nl- We compared results from this method 
with the geometric approach (see Ref. [13] for early work in this area; see Ref. [14] for some 
more recent work) used in our previous paper [6] , and they were found to agree for all models 
tested, as well as giving the same distributions of observables when tested on our landscape 
model. 

With regard to calculating /nl, compared to other methods (for instance Ref. [15, 16]), 
the Transport approach has the benefit of being computationally more efficient, as well 
providing a new means of understanding contributions to /nl by having explicit source 
terms. Equivalent to all other methods in the literature (including cosmological perturbation 
theory), it is simply an implementation of the separate universe assumption, but instead of 
evolving many perturbed trajectories, as is done in the popular 6N approach [15-17], one 
evolves probability distributions. 

What follows is largely a summary of the work done in Refs. [11, 12]. We focus on 
explicitly showing how the Transport formalism is implemented for a general d-field model 
of inflation and refer the reader to Refs. [11, 12] for the details. 

3.1 Moments of ( 

In calculating the statistical properties of the curvature perturbation we invoke the separate 
universe assumption and consider a collection of space-time volumes whose mutual scatter 
will ultimately determine the microwave background anisotropy on a given scale. Each space- 
time volume follows a slightly different trajectory in field space, whose position at a given 
time we label (/>*, the scatter of which is determined by the vacuum fluctuations at horizon 
exit. Here and in what follows the superscript "*" indicates that the quantity is evaluated on 
a spatially-flat hypersurface. If we know the distribution P{4>*) then, among other things, we 
can study the statistical properties of the deviation of these trajectories from their expectation 
value 6(j)* = (p* — where i indexes the components of the trajectory (f), namely the 
species of light scalar fields. The two-point correlations among the 5(j)i are expressed by the 
covariance matrix S(t), where 



The covariance matrix, third moment and centroid $ are all functions of time, but in our 
notation we will be suppressing the explicit time dependence. 

A consequence of the separate universe assumption [17] is that the curvature pertur- 
bation evaluated at some time t = tc is equivalent on large scales to the perturbation of 
the number of e- foldings N{tc,t^:,x) from an initial fiat hypersurface at i = t*, to a final 
uniform-density hypersurface at t = tc, 




(3.1) 



and the third moment is given by 




(3.2) 



C{tc, x) ~ 6N{tc, U,x) = N{tc, U,x)- N{tc, Q 



(3.3) 
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where 

N{tc,U) = Hdt. (3.4) 
Expanding 6N in terms of the initial field perturbations to second order, one obtains 

C(tc, x) = 6N{tc, t„x) = Ni5<P* + ^Ny{5(j)i6(j)j - {dc^iScPj)), (3.5) 

where repeated indices should be summed over, and N^i, N^ij represent first and second 
derivatives of the number of e- folds with respect to the fields 4>*. We remind the reader that 
it is necessary to subtract the correlation function in the second term. This is because one 
can interpret the covariance matrix as the contribution from disconnected diagrams which 
gives the vacuum energy. In Fourier space one only considers connected diagrams from the 
outset and thus the subtraction is already implicitly taken care of. 

Combining Eq. (3.1) and Eq. (3.2) with Eq. (3.5) we get expressions for the two- and 
three-point functions in terms of the moments of P{(f)*). The two-point function is 

{CO=NiNj^,j. (3.6) 

It is useful to decompose the three-point function as 

(CCC) = (CCC)i + (CCC)2, (3.7) 

where 

{CCOi = NiNjNkaijk, (3.8) 

and 

{CCO2 = ^N^NjNkm [^ik^jm + ^i^^jk] • (3.9) 

Eq. (3.8) is the intrinsic non-linearity among the fields, while Eq. (3.9) encodes the non- 
Gaussianity resulting from the gauge transformation to C; as one evolves from one fiat hy- 
persurface to another, turns in the trajectory will contribute to the non-Gaussianity. This, 
as well as any non-Gaussianity present at horizon exit, is what is encapsulated in Eq. (3.8). 
However, this super-horizon evolution also causes the hypersurface of constant density to 
change and so the gauge transformation from the flat hypersurface to the coinciding surface 
of constant density also contributes to the non-Gaussianity and this contribution is taken 
into account in Eq. (3.9). 

3.2 Derivatives of N 

From Eq. (3.8) and Eq. (3.9), it is clear that in order to calculate moments of the power 
spectrum we need a method for calculating derivatives of N. In general when using the 6N 
technique it is difficult or impossible to find an analytic expression for the derivatives of N. It 
is therefore necessary to run the background field equations many times from perturbatively 
different initial conditions, stopping at some value for H which is the same for all the runs. 
One then calculates the derivatives of N with respect to the initial conditions. In using the 
Transport equations, however, this process is replaced by solving a set of coupled ordinary 
differential equations. Instead of taking the surfaces "*" and "c" to be at horizon crossing and 
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time of evaluation^ respectively, instead the surfaces are taken to be infinitesimally separated 
and the transport equations evolve the field values at horizon crossing forward to the time 
of evaluation. The upshot of this is two-fold. As we will see, the use of ordinary differential 
equations to evolve the moments of the field perturbations allows us to see the source of 
super- horizon evolution and hence the various contributions to /nl- The second and more 
immediate benefit to our current discussion is that we can find a general expression for the 
derivatives of A^. To leading order in slow-roll, for a given species "z" , the number of e-folds 
N between the fiat hypersurface and a comoving hypersurface is given by 



N{tc,U) 



V 



no sum on i . 



(3.10) 



and so if the two surfaces are infinitesimally separated, then we can write 



dN 



V 



(3.11) 



* the method used in Refs. [16, 18] for sum-separable potentials is also now 



To handle 

applicable and we introduce the quantity 



+ 



V 



i+l 



(3.12) 



This enables us to write 



dC 



3 '^V^fc 



(3.13) 



which after some algebra gives the expression 



d(t)'i 



V 



(3.14) 



Hence we find 



V 



Vi 



VkV^ 



ky,k 



no sum on i 



(3.15) 



and 



AT - ^^^J j_ ^K^i _ 2VVikVkVj _ 2VVjkVkV,i ^ 2VViVkVkiViVj 



VkVk VkV 



ky,k 



(VkVk)'' 



iVkVk 



iVkVk? 



(3.16) 



where in Eq. (3.16) the limit c — )• * has been taken. 



''Time of evaluation is often taken to be the end of inflation but with regard to calculating observables, 
any time after isocurvature modes have decayed away will give the same result. A problem arises when 
isocurvature modes are still present at the end of inflation. In this case the power spectrum will continue 
to evolve and without a model of reheating this renders the model non-predictive. We will return to this in 
greater detail later on. 
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3.3 Transporting the moments 

Finally, we need a method for evolving the moments of the scalar perturbations. The proba- 
bility distribution P{4>*) is conserved and so, as described by the standard continuity equa- 
tion, the rate of change of P is given by the divergence of the current, 

where Ui = 4>i is the field velocity. The key achievement of Ref. [11] was to develop a method 
for extracting the evolution equations of the moments of P from the continuity equation. 
In Ref. [12] an alternative method was introduced, generalising to any number of fields on 
arbitrary slicing. We do not go into the techniques here; instead we just quote the resulting 
evolution equations for the centroid, variance and skew which collectively we refer to as the 
Transport equations, 

^'i = 'Pi + -^Ui^rnn'^mn + ••• , (3.18) 

^ij — fJ'i,m^mj ~l~ f^j,m^mi ~l~ ~^'^i,mn(^jmn ~l~ '^'^ j ,mn(^imn j (3.19) 

«'iife = Ui^mamjk + Ui^mn^jm^kn + (cyclic Z j fc) H . (3.20) 

The equation for the centroid Eq. (3.18) says that the mean field value evolves as the velocity 
of the fields but can be affected by evolution of the wings of the distribution. The evolu- 
tion equations for the variance and skew, as one might guess from the continuity equation 
Eq. (3.17), give evolution as the divergence of the field velocity but now also with source 
terms coming from the other moments. 

3.4 Cross-sections of the bundle 

As will be discussed in more detail in due course, an important consideration in our analysis 
will be whether or not evolution of observables is still taking place at the time of evaluation. 
Evolution stops when the trajectory becomes effectively single field [19]. This is to say the 
trajectory has reduced to a caustic [20] , so for this reason we would like a description for the 
evolution of the cross-section of the perturbed trajectories. Such a description has recently 
been developed in Ref. [21], to which we refer the reader for more detailed discussion. For 
simplicity we only describe the broad concept here and quote results that will be needed in 
future discussion. 

Cross-sections within the bundle are focused, sheared and rotated by the fiow. These 
distortions can be characterised by the evolution of connecting vectors describing the dis- 
placement between nearby trajectories in the bundle. If 6xi is an infinitesimal connecting 
vector, then assuming Uij is sufficiently smooth, 6xi is transported as 

6xi = 5xj—-^ (3.21) 

It follows that changes in the cross-section of the bundle can be determined in terms of the 
expansion tensor Ujj. We can decompose this in terms of a dilation 6 = tr Uij, a traceless 
symmetric shear fjjj, and a traceless antisymmetric twist Uij, 

Q 

Uij = -6ij + aij + LOij (3.22) 
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Dilation describes a rigid rescaling of 6xi by 1 + 6', representing a global tendency of the 
trajectories to focus or defocus. The shear encapsulates the tendency for some trajectories 
to flow faster than others while conserving the cross-sectional area of the bundle. The twist 
represents a rotation of the bundle with preserved volume, such as the tendency of trajectories 
to braid. The dilation, shear and twist act as sources for one another and so one expects a 
bundle will typically exhibit all of these behaviours at some point. 

We refer the reader to Ref. [21] for a more formal description of this formalism and its 
applications, but for the purposes of this paper all we need is the result that the focussing 
of the bundle is given by 

^1 



Q{H,Ho) = exp 



d 



H 



e{h)dh 

Ho 



(3.23) 



4 Findings 



Having set up our models and the machinery necessary to compute the observables, we now 
proceed to our results. The principal variables of interest to vary are the number of fields d 
and the horizontal mass scale mh. Large values of the latter correspond to relatively smooth 
potentials, and small values to heavily featured potentials. We refer to individual realizations 
giving sufficient inflation as 'verses'. 

We discuss our results in the following sequence: 

1. Dynamical properties of trajectories. 

2. Perturbation evolution along individual trajectories. 

3. Distribution of observables over ensembles of trajectories. 
4.1 Trajectory dynamics 

Minimal d dependence Fig. 3 summarises the qualitative behaviour found during our 
exploration of the properties of multi-field trajectories. We see that A(f), the length of the 
trajectory in field space, and i?^, the percentage increase in Acj) due to turns in the trajectory, 
defined as 

= im— — ^ — ^— ^ (4.1) 



show only a mild sensitivity to the number of fields. In fact, for all observables and infla- 
tionary parameters we looked at, the sensitivity to changing the number of fields was small 
over the range d = 2 to d = 6 compared to the spread of results. This was true even for 
the relatively predictive large-field case of mh = 15Mpi. Least sensitive of all, we found 
no change whatsoever in the distribution of slow-roll parameters at horizon crossing with 
e* = 0.002 ± 0.002 and r]* = 0.0012 ± 0.0004. 

e-fold distributions One instance where we did find sensitivity to d was in the e-fold 
distributions, where we saw a decrease in the proportion of trajectories with more than 60 
e-folds from 0.08 for two fields down to 0.01 for six fields. Fig. 3 hints at a tendency for 
trajectories to become shorter and more curved as the number of fields increases. So given 
that the trajectory will seek the route of steepest descent, it appears that for our model, 
increasing the number of fields increases the chance of the trajectory encountering a slope 
sufficiently steep to kill inflation. 
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Figure 3. How the mean A(/) and percentage increase in due to bending of the trajectory changes 
with the number of fields, for toji — 15Mpi. The bars show how the standard deviation of these 
quantities changes over this range. 



In Refs. [1, 3] slightly different parameterisations of the inflationary region of the po- 
tential were used to show the probability of obtaining a given number of e-folds of inflation 
was 

PiN) oc ^, (4.2) 

where a = 4 and a = 3 were found respectively. Bearing Eq. (2.3) in mind, it seems 
reasonable that Eq. (4.2) might apply more generally and indeed we find it to be a good fit 
to our e-fold distributions, with a increasing as the potential becomes more featured from 
a = 2 through to q = 5 over the range of mh we investigated. To illustrate the implication of 
this, consider the popular idea that inflation was preceded by a tunnelling event. Using the 
(perhaps somewhat arbitrary) values of Ref. [3], we place an anthropic lower bound on the 
number of e-folds at 59.5 coming from structure formation and an observational lower bound 
at = 62 on the curvature from tunnelling. Then for a = 2, the probability of the model 
achieving sufficient inflation to be in agreement with observation is roughly 92%, while for 
a = 5 it is more like 81%. However, remember we are working in the range of A(p which is not 
best motivated theoretically. While not accessible with the techniques used here, our results 
lead us to expect that for small-field models, a should be larger. This has the potential to 
cause tension with observation, as by a = 18 the chance of finding ourselves in the observed 
universe falls to 49%, i.e a typical observer would expect to see evidence of curvature. 

4.2 Perturbation evolution: P^^, /nl and the adiabatic limit 

Having seen that our model is insensitive to the number of light fields, for simplicity we only 
give results for two-field potentials in the remaining sections of this paper, focussing mainly 
on the dependence on the feature scale of the potential. But we would like to emphasise that 
the results hold more generally. 

4.2.1 and 9 

The key difference between single-field and multi-field inflation is that the latter admits 
evolution of the power spectrum on super- horizon scales. This means that in order to make a 
prediction from multi-field inflation one needs to know the full evolutionary history up until 
the model becomes effectively single field, i.e. the adiabatic limit is reached [19, 20]. Once 
this happens the power spectrum stops evolving and one can evaluate observable quantities 
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Figure 4. Example plots showing the super-horizon evolution of the power spectrum for modes 
exiting 55 e-folds before the end of inflation, mh = 3Mpi corresponds to a highly-featured potential, 
and in all such cases evolution stops long before the end of inflation, mh = 15Afpi is a comparatively 
smooth landscape, and then a significant proportion of verses are still evolving at the end of inflation. 



at a subsequent time of one's pleasing. The problem is that there is no guarantee that such an 
adiabatic limit will be reached before the end of inflation, and if this is not the case making 
a prediction requires knowledge of reheating and so forth. 

As can be seen in Figs. 4, 5 and 6, the ability to reach the adiabatic limit, where the 
perturbation on a given scale becomes constant, is strongly dependent on how featured the 
landscape is. For the more featured landscapes such as mh = 3Mpi, we found the adiabatic 
limit was reached in all cases (even the trajectories disappearing off the top of the plots), 
while for the smoother landscapes like mh = 15Mpi the proportion of trajectories achieving 
this clearly drops significantly. 

There is a very intuitive reason for why this should be the case. Rewriting Eq. (3.23) 
in terms of e-folds N we find 



f-N 

e{N,N^) = exp ' ' 



1 

- / {3e - 2f] + ti Mij)dN 
where Mij is the Hessian of In V and f/ is the generalised slow-roll parameter 



(4.3) 



We therefore see that in a valley, strong focussing will occur, while on a ridge or a hilltop the 
bundle will dilate. With this picture it is quite easy to see why we should expect evolution 
as seen in Fig. 4. The treacherous landscape of rrih = 2>M-p\ typically gives exactly the 
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Figure 5. Example plots showing the super-horizon evolution of /nl for modes exiting 55 e- folds 
before the end of inflation. Again we see that for the more featured landscape evolution stops early 
on, while for the smoother example, evolution often continues to the end of inflation. 

100 

so 

60 
40 

20 

. nllnin aan - .n i .- . i ■ — . l @ 

02 0,4 0,6 0.8 1.0 

Figure 6. Q at the end of inflation for mh = 3Mpi (blue), mh = 9 Afpi (yellow) and Wh — 15 Afpi 
(red). All trajectories essentially reach a caustic in the most featured example, less for mh — 9A/pi 
and least for the smoothest landscape mh = 15Mpi. 

conditions required for a very strong focusing, while in contrast the comparatively mild, 
undulating meadows of mh = 15Mpi give very little incentive for trajectories to focus to a 
caustic. 

4.2.2 /nl is always small 

Much of this kind of discussion carries over to understanding the results of Fig. 5. First 
and foremost it should be noted there was not a single example of a trajectory that gave 
sufficient non-Gaussianity to be detected by any future planned experiment. Methods to get 
around this disappointingly generic feature of multi-field inflation were recently addressed in 
Ref. [20] and the special case of sum-separable potentials was also discussed in Ref. [22]. In 
the case of sum-separable hilltop potentials which reach an adiabatic limit during inflation, 
what is known as the horizon-crossing approximation [23] gives a good estimate of the final 
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Figure 7. Plots of Verse 113147 (left) and Verse 253911 (right) showing very distinctive /nl evolu- 
tion. In the plots on the top row, /nl is blue, the intrinsic component, /nli, is red and the gauge 
contribution, /nl2) is yellow. 



value for /nl- 



/nl 



'6 14 



(4.5) 



where in this instance (j) represents one or at most a few fields where Ni is large. From this 
we see that provided there are enough fields to keep jy small, with the right initial conditions, 
a ridge can give rise to a large /nl- 

What we find in our analysis is that the problem of obtaining a large /nl is made 
particularly acute by the need to obtain sufficient inflation. When we have a very smooth 
landscape, sufficient inflation is easily achieved but the lack of features means there is nothing 
to give rise to a large /nl- On the other hand, when the potential is very featured, it is 
difficult to start close enough to a ridge to get a large /nl without falling off it, thereby 
killing inflation. In some models, such as axion N-flation [24], a sufficiently large number of 
fields can make it possible to overcome this problem [25] , as the large damping term makes it 
easier to be close to a ridge without falling off too soon. For our model, while we do not have 
the computational power to explore this possibility, with enough fields we would expect to 
see some examples with a large /nLi but they would constitute only a very small proportion. 
This is because if many fields have a large N^i, their contributions to /nl will, through a 
manifestation of the central limit theorem, cause /nl to be vanishing in the limit of many 
contributing fields. Thus on average we would always expect /nl to be small. 
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4.2.3 Trends in /nl evolution 

While we found no examples of large /nl we did find a very diverse range of behaviour. 
This diversity is indicative of why detection would be such a powerful constraint on models. 
That said we did find some common trends. Fig. 7 shows two examples of evolution of /nl 
together with the corresponding evolution in the width of the bundle. We chose these two 
examples, Verse 113147 and Verse 253911, in particular as each shows characteristics that 
were common to most trajectories, but also demonstrate that counter examples were found. 

The gauge contribution determines the peak As particularly well demonstrated by 
Verse 113147, in all examples we found /nl2 determined any peaks in /nl- Typically it was 
the case that the intrinsic non-Gaussianity played a highly subdominant role in the feature, 
but we did find the exception that is Verse 253911 which clearly received an important 
contribution from the intrinsic part. In Ref. [11] it was noted that for the double quadratic 
potential and quadratic exponential potential this behaviour was present. Here we show this 
characteristic applies much more generally. 

/nl grows when the bundle dilates Features in /nl occur whenever Q grows. This seems 
very reasonable since at this point the perturbed trajectories will be exploring different parts 
of the potential. Typically the peak in /nl occurred very close to the time of the peak in Q 
but again. Verse 253911 shows this need not be the case precisely. We see features in /nl2 
are intimately related to features in 9. We attribute this to common terms involving V^ij. 

Asymptotic behaviour of /nl is not straightforward A result that continues to elude 
us is a simple way of understanding what the final value of /nl in the adiabatic limit will 
be. For sum-separable hilltop potentials, the horizon-crossing approximation works well, 
but for more general potentials there is no equivalent. As the examples in Fig. 7 show, the 
asymptotic value can be reached in dramatically different manners. Worse still is the fact 
that in the adiabatic limit the intrinsic and gauge transformations need not settle to constant 
values. This indicates that a different set of parameters should be considered if we are to 
make progress with this question. 

4.3 Interlude: The Lyth bound 

Before moving on to discuss distributions of observables, we would like to take a brief moment 
to discuss the relation between field trajectories and the tensor-to-scalar ratio, as it will be 
helpful to bear in mind in the subsequent discussion. 

Taking A'^cmb to be the number of e-folds between when fluctuations on CMB scales 
left the horizon and the end of inflation, we can obtain a d-field version of the Lyth bound 
by writing 



where we are integrating along the field trajectory. If we assume e is either constant or 
increasing over this period, then we have 




(4.6) 



2e < 



(4.7) 



A^CMB 



For single field inflation r = 16e but when there are more fields the curvature perturbation 
evolves on super-horizon scales, suppressing r and so r < 16e [6]. We therefore see that the 



Lyth bound remains essentially the same for multi-field models as in the single-field case 



Planck hopes to measure the tensor-to-scalar ratio with an accuracy of a few hundredths, 
hence has discovery potential if it is of order 0.1 or so. Comparing the Lyth bound with 
the discussion in section 2.3 we therefore see that a detection would exclude all small- and 
medium-field models if only one field is admitted. However as previously discussed, multi- 
field models are strongly motivated by fundamental theory. If we consider the extreme case 
of sum-separable potentials then the discussion of Section 2.3 requires each A0j ^ Mpi but 
there is no restriction on the number of fields contributing during infiation; this was for 
instance the motivation of the N-flation proposal [24]. Therefore, if we are to stay in the field 
theory favoured regime of small- field models, a detection of r would place a lower bound on 
the number of fields! Rewriting Eq. (4.8) in a more suggestive form we have 



Fig. 8 shows our findings for n and r at the end of infiation, also summarised in Table 2. As 
the landscape becomes more featured the viable trajectories become shorter and increasingly 
bendy. The Lyth bound tells us the distance travelled in field space places an upper bound on 
the tensor-to-scalar ratio. Furthermore, bends in the trajectory cause super-horizon evolution 
of the curvature power spectrum, while the tensor power spectrum is conserved, and so as 
Fig. 8 shows, we see an increasingly strong suppression in r as we move to lower mh. 

As summarised in Fig. 6, a less featured potential reduces the chance of trajectories 
reaching their adiabatic limit. It is noteworthy that, despite this, the plots of the n-r plane 
for mh = 15Mpi and mh = 18Mpi show remarkable consistency for n and r at the end of 
infiation. This might lead one to think that the super-horizon evolution is having negligible 
effect, but if we take the example of mh = 15 Mpi, as Table 2 shows, the mean increase in the 
field trajectory from turning is only 4%, yet if we assumed a single- field approximation was 
valid we would obtain n = 0.98 ±0.01 which is significantly different from the actual result 
of 0.95 lb 0.01. If nothing else, these results show one should be exceedingly careful when 
making single-field approximations. Fig. 9 compares distributions for the spectral index with 
those obtained using a single-field approximation. 

The central concern regarding the possible existence of a landscape is whether or not 
such a model can be tested. As we have mentioned, a key challenge is the measure problem, 
another being our very limited understanding of fundamental theory. Once these problems are 
better understood though, we will still be left with distributions for observable quantities. 
No matter how well developed our understanding, it seems reasonable to assume that at 
some level our ability to make predictions will be fundamentally limited by the details of the 
theory. Our toy landscape illustrates this in a very explicit way. In a sense mh gives a way 
of quantifying the complexity of each landscape. For our model, we see the spread of results 
for the spectral index dramatically increases as we move to more featured landscapes, while 
for the tensor-to-scalar ratio the spreading is considerably less dramatic due to a suppression 
coming from the inevitable decrease in the length of the field trajectory A(j). As we have seen 
/nl by contrast remains consistently small. 







4.4 Distribution of observables: n and r 
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Figure 8. n-r plots for a range of TOh, beginning witli the more featured potentials. The curve shows 
the WMAP7+aU 95% confidence limit [26]. 



mh(Mpi) 




B4%) 


n 


r, 95% conf. 


/nl 


3 


4.3 ± 1.7 


17± 18 


1.03 ±0.15 


0.036, r < 0.1981 


-0.008 ±0.078 


6 


7.6 ± 1.8 


16± 18 


0.93 ±0.07 


0.036, r < 0.172 


0.029 ± 0.034 


9 


8.4 ± 1.8 


10± 13 


0.93 ±0.03 


0.023, r < 0.081 


0.028 ±0.022 


12 


8.6 ± 1.4 


5±7 


0.94 ±0.02 


0.024, r < 0.066 


0.022 ±0.010 


15 


9.0 ± 1.3 


4±7 


0.95 ±0.01 


0.032, r < 0.074 


0.018 ±0.007 


18 


9.2 ± 1.2 


4±9 


0.96 ±0.02 


0.039, r < 0.079 


0.016 ±0.009 



Table 2. Table of the mean distance in field space travelled in the last 55 e-fold of inflation for a 
given rrih, B^, the mean percentage increase coming from bends in the trajectory, and corresponding 
results for observables. 

5 Conclusions 

We explored inflationary dynamics in randomly-generated potentials as well as the con- 
sequences for super-horizon evolution of perturbations. We found this exploration to be 
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Figure 9. Example plots comparing spectral index as obtained using single-field approximation (red) 
and as calculated taking evolution up to the end of inflation into account (blue). 



interesting primarily on two fronts. 

First, by exploring a very large number of inflationary trajectories, we encountered a 
wide range of super- horizon evolution behaviour for P(^(^ and /nl- The benefit of this was 
that it became easy to see what characteristics are generic and which are not. We found 
that peaks in /nl tend to be determined by the gauge contribution but behaviour was rather 
more broad in the adiabatic limit, showing few trends. Understanding of non-Gaussianity 
is still in rapid development and so exploration of this kind can be very helpful in gaining 
insight in how to progress towards something more concrete such as Ref. [21]. 

We also found that keeping track of the easy-to-compute bundle width was extremely 
informative. By following how G changes along the trajectory, we were generally able to 
understand what qualities of the potential gave rise to super-horizon evolution. In particular, 
we found that peaks in /nl occur during regions of the potential that give rise to a dilation 
of the bundle. However, a more quantitative description awaits future development. We 
emphasise that in order to make predictions in any multi-field model, one needs to perform 
an equivalent analysis to ensure no evolution is taking place at the time of evaluation. We 
found that as the mean length of the field trajectory in field space increased, the chances of 
reaching an adiabatic limit drastically decreased, rendering the larger field models essentially 
non-predictive without a model of reheating. 

Second, we looked at how varying the scale of features and the number of light fields 
affected the ensembles produced for a given parameter. We found that landscapes where the 
mean length of the field trajectory was large typically gave results consistent with current 
observational data (despite not necessarily reaching their adiabatic limit). However for more 
featured landscapes where the mean field trajectory was smaller, the spread in the spectral 
index increased significantly. The spread in the tensor-to-scalar ratio did not increase so 
dramatically. This can be understood in terms of the Lyth bound which places an upper 
bound on the tensor-to-scalar ratio according to the length of the field trajectory. We found 
varying the number of fields between 2 and 6 to have negligible effect on the distributions 
for observables. Amongst all trajectories we found no examples of detectably large non- 
Gaussianity. In absence of motivation for why we would be an atypical observer, this result 
is sufficiently strong to conclude that an observation of local-type non-Gaussianity by Planck 
would rule out models of this kind. 
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